Theory of High Tc Ferromagnetism in SvBq family: 
A case of Doped Spin-1 Mott insulator in a Valence Bond Solid Phase 
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Doped divalent hexaborides such as Sri-^La^B^ exhibit high Tc ferromagnetism. We isolate a 
degenerate pair of 2p-orbitals of boron with two valence electrons, invoke electron correlation and 
Hund coupling, to suggest that the undoped state is better viewed as a spin-1 Mott insulator; it 
is predicted to be a type of 3d Haldane gap phase with a spin gap ~ 0.1 eV , much smaller than 
the charge gap of > 1.0 eV seen in ARPES. The experimentally seen high Tc 'ferromagnetism' is 
argued to be a complex magnetic order in disguise - either a canted 6-sublattice AFM (« 120°) order 
or its quantum melted version, a chiral spin liquid state, arising from a type of double exchange 
mechanism. 



The observation ||l| of high 'ferromagnetism' in 
lightly doped SrBg family is a great surprise in con- 
densed matter physics in recent times; neither Sr nor 
B are known to participate in magnetism. Another sur- 
prise that followed was a high Tc superconductivity in 
MgB2 , a diboride. Electron deficient B is known to 
form molecules and solids with varying ligancy, from the 
stable doubly charged octahedral {BqHq)'^~ molecule to 
a metallic phase of boron, that has B12 icosahedral clus- 
ter as a basic building block. How a 2p atom like boron in 
solid state manages to achieve high Tc superconductivity 
as well as high Tc ferromagnetism is a fascinating ques- 
tion. 

There is sustained effort to understand ferromagnetism 
in the SvBq family. Some of the new experimental results 
are intriguing and can not be explained in a satis- 
factory fashion by the existing theories including 
the popular excitonic instability [^-13| scenario. Notable 
among them are the recent ARPES and x-ray emission 
studies 0, that show a gap > 1 eV at the Fermi level 
for pure SrBg . A gap of ^ 0.8 eV is also predicted by 
an electronic structure calculation ||l^ using GW method 
that attempts to incorporate correlation effects in earlier 
approaches ]l5|-Jl7[|. 

In the present letter we propose that it is advanta- 
geous and perhaps also correct to view the StBq family 
as a spin-1 Mott insulator, in view of strong coulomb re- 
pulsions and Hund coupling among two electrons in two 
degenerate valence orbitals of boron. However, at the 
outset we should point out that the difference between 
a band and Mott insulator, in an even electron number 
(per unit cell) insulator (with no magnetic order) like 
ours, is quantitative in the sense it is largely determined 
by whether Et << Ec or not. Here Et is the spin gap 
or the lowest magnetic triplet exciton energy and Ec is 
the charge gap(as seen by ARPES for example). 

Our present proposal is consistent with all the known 
experimental results and further certain novel predic- 
tions, which can be experimentally tested, follow in a nat- 
ural fashion, i) SrBg is a spin-1 Mott insulator, with two 



electrons in a doubly degenerate 2p-orbitals providing a 
spin-1 moment, ii) Antiferromagnetic coupling among 
the localized spin-1 moments, arising from kinetic ex- 
change leads to a type of 3d Haldane gap or Majumdar- 
Ghosh phase with unbroken lattice symmetry{&gvLTe 1). 
The singlet valence bonds are 'frozen' at the B — B bond 
bridging two neighboring Bq octahedra. This phase has 
a small spin gap Et ~ 0.1 eV (this is not incompatible 
with a diamagnetic behavior of the insulating SrB^ )■ 
iii) Doping liberates, through a form of double exchange 
mechanism, a canted 6-sublattice AFM (« 120°) order 
or a chiral spin liquid state, {Sia ■ {Sip x 8^-^)) 7^ both 
with a small ferromagnetic moment. Our predictions 
could be tested by neutron scattering or other low en- 
ergy probes. Two of our robust predictions are a small 
spin gap, Et ^ 0.1 eV and at the least a well developed 
120^^ spin correlations (figure 1) inside a Bq octahedron. 

The following two experimental facts about SrB^ , 
when taken together are striking and gives an impor- 
tant clue for our model building, i) As revealed by recent 
ARPES experiments the parent insulator has a charge 
gap > 1 eV ii) The insulating paramagnetic ground state 
is very fragile (unlike a band insulator) and a small dop- 
ing in S'ro.ggs-Lao.oosSei that adds half a percent of car- 
rier per La atom, produces a ferromagnetic phase (with 
a small moment w O.OT/is per La atom) with a large 
Tc ~ 600 - 900X. 

Some key quantum chemical information about the 
SrBe family that we will use in building our theory are: i) 
each B atom has two valence electrons and two degener- 
ate valence p-orbitals, ii) the intra atomic Hund coupling 
is w l.deV, iii) the unscreened Hubbard U for the 2p or- 
bitals of B can be as large as 8 — 10 eV. iv) The nearest 
neighbor inter octahedral B — B distance is smaller than 
that inside an octahedron by about 5 percent. 

We first briefly present the conventional electronic 
structure ||l^ description of Bq cluster in SrB^ . Bq 
octahedra are covalently bonded to form a simple cubic 
lattice. There is nearly complete charge transfer from 
Sr to ^6 clusteriSr^+^e^". I gnoring the core orbitals of 



1 



B as well as Sr'^'^ we are left with one 2s and three 2p 
orbitals per B atom. There are 20 electrons in these 24 
B orbitals in every octahedron. The 24 orbitals can be 
separated into two sets: i) sp^ hybrids (two per B atom) 
that are along the body diagonal of the Bq octahedra 
and ii) two p-orbitals per B atom that are tangential to 
the octahedra (inset in figure 1). An s-like combination 
of the six sp^ hybrids that are pointing at the octahe- 
dral center form a very strong six center bond with a 
binding energy of about 15 eV. This bonding state takes 
two electrons and is primarily responsible for the stabil- 
ity of the octahedra. The sp^ orbital pointing radially 
outwards strongly hybridize with the corresponding or- 
bital of neighboring octahedra and results in a stable co- 
valently bonded cubic network of Bq octahedra. These 
bonding states take away 6 electrons per Bq octahedra. 
We are now left with 12 electrons and 12 p- orbitals per 
octahedron, i.e., two tangential 2p orbitals and two elec- 
trons per B atom. 

In the first band structure calculation for 
SrB^ family, Longuet-Higgins and Roberts, using the 12 
p-orbitals, form a set of four triply degenerate molecular 
orbitals - tiu,t2g,tig and t2g. (For simplicity we will ig- 
nore the hybridization of tin with symmetry adapted sp^ 
orbitals that leads to t[^ and t"iu). These Bq cluster 
orbitals overlap to produce 6 bonding and 6 anti bond- 
ing Bloch bands. The 6 bonding bands are completely 
filled by 12 electrons to produce a band insulator. Later 
workers p^ , p^ emphasized the mixing of anion d-states 
and showed that within band theory it reduces the band 
gap to nearly zero value at the X-points in k-space. 

In the absence of any electron-electron interaction the 
energy difference between the top most tig and and bot- 
tom most tiu molecular orbitals is ~ 10 eV. This is a 
measure of the kinetic energy of delocalization per elec- 
tron in the slater determinant state within a Bq cluster. 
This energy is comparable to the Hubbard U ~ 8 — 10 eV. 
In addition we have a ferromagnetic Hund coupling be- 
tween two p-orbitals of a i? atom, Jh ~ 1.5 eV. As 
the kinetic energy of delocalization within the cluster is 
comparable to the energy increase from Hubbard U and 
Hund coupling Jh, the low lying cluster eigen states are 
strongly perturbed by many body effects, to the extent 
the simple Bq molecular orbitals and their slater deter- 
minant states mostly loose their relevance. 

Thus a natural starting point to understand the low 
energy physics of the cluster is the Mott localization of 
two electrons and form a spin-1 boron moment. Our ap- 
proach is similar to restricting oneself to the valence band 
basis in the case of carbon pi: bonded system such as ben- 
zene (with very similar quantum chemical parameters), 
where it is known to work very well; ours is a general- 
ization of valence bond basis approach to the case of two 
p-orbitals per site. Further the inter Bq cluster hoping 
does not modify this localized picture significantly, as the 
inter molecular orbital hopping matrix element between 



neighboring clusters is small '--^ 0.25 eV (that leads to a 
band width of 12 x 0.25 w 3 eV, as seen in band structure 
results). 

We have formalized the above by starting from a dou- 
bly degenerate 3 dimensional Hubbard model containing 
an average of one electron per orbital: 

— «/_ff ^ Siai • SiQ2 (1) 

Here t'^^'^'^ represents the nearest neighbor hopping in- 
tegrals; i denotes the centers of octahedra that form a 
simple cubic lattice, a — 1,..6 denotes a B site within 
an octahedron and pi — 1,2 denotes the two degenerate 
p-orbitals. The operator s^q^ is a spin half operator of 
an electron in the /ith p-orbital of a'th B site in the i-th 
octahedron. 

Our Hubbard model contains three types of near- 
est neighbor hopping matrix elements t-j^i, tT^2, and iTro- 
shown in the inset in figure 1; their values range over 
1 eV to 1.5 eV . As mentioned earlier U ~ 8 eV and 
Jh 1.5 eV. 

To understand the low energy spin dynamics we per- 
form a kinetic exchange perturbation theory and get a 3 
dimensional spin-1 Hamiltonian for our Mott insulator: 

Hs — Jl ^ Si^a ■ Si,/3 + J2 ^ Si^a ■ Sj,l3 (2) 

i,(a/3) 

The first term represents the Bq cluster spin Hamil- 
tonian and second the octahedral bridge spin coupling. 
Here S^q is the spin-1 operator of the a'th B atom in 
the i-th Bq octahedron. Ji and J2 are intra cluster 
and inter cluster kinetic exchange integrals. A direct 
kinetic exchange perturbation theory gives a large an- 
tiferromagnetic coupling and it needs to be corrected by 
subtracting a direct or 'potential exchange' (arising from 
non-orthogonality of the nearest neighbor orbitals) that 
favors ferromagnetic alignment. We have estimated the 
potential exchange by looking at some experimental re- 
sults and quantum chemical calculations for the related 
carbon systems. After making this subtraction we esti- 
mate Jl, J2 '-^ 0.1 to 0.3 eV and Ji < J2. In view of the 
approximate nature of our estimates, we will keep Ji and 
J2 as parameters to be determined experimentally. 

Now we derive the phase diagram for our spin Hamil- 
tonian as a function of its only parameter, the dimen- 
sionless ratio (figure 2). It is convenient to rewrite 
the Hamiltonian (equation 2) as 

~ ~^ (Si,Q + ^i,p + 

+ Y II (S,,a + S,,^)2-7V(12Ji-t-6J2), (3) 
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FIG. 1. Schematic phase diagram of our Boron spin-1 
Hamiltonian. Doubles lines denote 'frozen' valence bonds at 
the octahedral bridges. Inset shows the 'tangential' boron 2p 
valence orbitals and 3 types of hopping matrix elements. 

where N is the total number of octahedra. The first 
term of equation (3) represents the sum over triangles 
of three spins forming the eight faces of a Bq octahe- 
dron. Each term of equation 3 is a positive operator 
with eigen values > 0. It follows immediately that the 
classical ground state of the above Hamiltonian is a 120" 
six-sublattice antiferromagnet for any Ji, J2 > 0. The 
spin pattern within an octahedron is shown in figure 2. 
Spins pairs at opposite corners of an octahedron are par- 
allel. Three such diagonal pairs within an octahedron are 
coplanar and at an angle of 120", among themselves. This 
defines a spin chirality (vorticity) ±1 for an octahedron. 
The cubic lattice of octahedra form two cubic sublattices 
with corresponding spins exactly anti-parallel, making it 
a 6-sublattice planar antiferromagnet (figure 1). The net 
result is that the bridge spin coupling is not frustrated, 
but the twelve nearest neighbor coupling within an oc- 
tahedron are frustrated. In spite of the spin reversals 
the chirality has the same sign in all the octahedra; i.e., 
we have a ferromagnetic chiral order. Thus our ground 
state exhibits a discrete two fold chirality (Ising like) de- 
generacy in addition to the global SU(2) spin rotational 
degeneracy. 

For ~ 1 the classical ground state exhibited above 
is stable and survives spin wave fluctuations. However, 
when >> 1 or « 1 quantum fluctuations destabi- 
lize the classical ground state and we get singlet ground 
states with a finite spin gap. 

When ^ << 1 each Bq cluster has an unique singlet 
ground state. There is a finite gap ~ Ji for spin-1 exci- 
tation. The couplings J2 between neighboring octahedra 
through the bridge spins reduce this gap by virtual exci- 
tations. We have estimated by perturbation theory that 
the gap survives until ^ We call this phase as the 
octahedral singlet sohd (OSS). 



When 4^ >> 1 the AFM order is destabilized; the 
octahedral bridge pairs become non-degenerate singlets 
with a gap of ~ J2 for spin-1 excitations. Spin couplings 
within Bq cluster reduce this gap by virtual excitations 
and the gap vanishes when we decrease ^ to & value ~ 2. 
This valence bond solid phase is a 3 dimensional realiza- 
tion of spin-1 Majumdar- Ghosh phase, that also retains 
the lattice symmetry. Further, this phase is also con- 
tinuously connected to an AKLT type of 3d Haldane 
gap phase, with quantum fiuctuating effective spin-3 mo- 
ments at Bq clusters 

Where is SrBg family in the above phase diagram ? 
It is unlikely that undoped SrBQ family has long range 
AFM order. The shortness of the octahedral bridge B—B 
distance by about 5 percent, compared to nearest neigh- 
bor B — B distance within an octahedra seen in all the 
three hexaborides, SrBe , CoBq , BqBq could also arise, 
apart from some quantum chemical reasons, from the for- 
mation of a 'frozen' singlet (valence bond) between the 
B spin-1 moments of the octahedral bridge pair. So it is 
very likely that the strong quantum fluctuations of the 
spin-1 system and the fact that > 1 in the SrBQ family 
is keeping it in a valence bond solid phase. 

Now we discuss how very small doping leads to a fer- 
romagnetic order with a large curie temperature. The 
analysis of the doped situation starting from the Hub- 
bard model (equation 3) is hard compared to the un- 
doped case. However, our preliminary analysis indicates 
very rich possibilities, all arising from the Mott insula- 
tor parentage and the 'vicinity' to an antiferromagnetic 
order. 

We believe that the observed small moment gives us 
an important clue as to the origin of magnetism. A nat- 
ural way in which a small moment can arise is by canting 
of an existing antiferromagnetic order. Our valence bond 
solid phase, being close to the 6-sublattice AFM phase in 
the phase diagram, should have short range antiferromag- 
netic correlations. Further if doped carrier delocalization 
decreases the effective ^ of our boron spin Hamiltonian, 
we my be pushed into the 6-sublattice AFM phase: we 
argue below that this is likely to happen. 
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FIG. 2. Two Possible Schematic Phase diagrams in the x-T 
plane. The experimentally observed 'ferromagnetic' phase is 
suggested to be a canted 6-sublattice antiferromagnetic or a 
chiral spin liquid phase 
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In a double exchange process, an added electron 
prefers higher total spins such as 2 or 1 between neigh- 
boring B spin-1 pairs, in order to satisfy Hund's rule 
during delocalization; singlet pairs are bottle necks for 
doped carrier delocalization. In the VBS phase such sin- 
glet amplitudes among neighboring spin pairs within an 
octahedron are relatively low: total spin values 2, 1 and 
of a pair occur with probabilities ~ |,| and i; on 
the other hand the singlet amplitude is larger '--^ 1 at an 
octahedral bridge pair. Carrier delocalization thus will 
selectively project out large local singlets from the bridge 
pairs; i.e., it effectively adds a local term, a singlet pro- 
jector, xtT!-2iSia + ^jf3)'^ , to the singlet dominated octahe- 
dral bridge pair. Thus Ji and J2 of our spin Hamiltonian 
(equation 2) get modified in a state dependent fashion to 
effective J2{x) « J2 — a;t7r2 and Ji{x) « Ji. 

The above results in a linear decrease in the ratio with 
doping x: ~ 77 ~ ^TT ' "^^^^ doping may allow us 

to enter (in figure 1) the 6-sublattice AFM phase (with 
a small canting induced ferromagnetism, to be discussed 
below) or a quantum melted version with long range chi- 
ral order. The scale of maximum Tc in our picture is 
roughly the maximum scale of of the insulating part of 
the phase diagram, determined by Ji J2: ksTc^ Ji, J2, 
which can be easily as large as 900 K. 

Once we establish a well developed short range or long 
range AFM order canting is easily obtained in the dou- 
ble exchange mechanism, as explained by de Genncs in 
1960 in the context of doped manganites. This is ex- 
plained by minimizing the sum of the exchange energy 
of an octahedral bridge pair spins, for example, and the 
Anderson-Hasegawa double exchange term : 



E fa J2 cos 9 ~tTj2 X cos ■ 



(4) 



Here 



is the angle between the two spins. The energy 
" = i.e., e « 180° - 

rather than at 



minimum occurs at cos 2 — 4j , i-^-, ^ ~ 4j 
180°. This leads to a small moment 



of w 0.9 /Is per formula unit, close to what is seen ex- 
perimentally. 

The preexisting planar chirality (vorticity) in our 6- 
sublattice AFM order, elaborated earlier, gives us a 
novel possibility of a finite non-planar chirality order p^ ] 
{Sia ■ (Sj^ X Si-y)) 7^ 0, through quantum fluctuations 
arising from dopant carrier delocalization. How cant- 
ing, a net magnetic moment and a non-planar chirality 
arise through quantum fluctuations is nicely illustrated 
in the following example. Consider three spin-half mo- 
ments coupled antiferromagnetically: H = J{Si ■ S2 -I- 
S2 ■ S3 -I- S3 • Si) has a 120° classical planar ground state 
{Sz = 0) with a two fold planar chirality degeneracy. 
However, the exact ground states have non-planar chiral- 
ity (Si • (S2 X S3))g = ±2V3 and a net spin =i moment, 
both arising from quantum fluctuation induced canting. 
One should not also rule out the possibility of nematic 



order such as (Sq, x Sp)q 7^ in view of the complexity 
in our system. 

Thus canting, weak moment ferromagnetism and de- 
velopment of non-zero chiral order parameter (S^q, • (S^^ x 
Si^)) ^ are all tied together. As the energy scales and 
the origin of the chirality stiffness and the 6-sublattice 
AFM stiffness are different the career delocalization can 
also quantum melt the long range AFM order but leave 
the chirality order intact. In this case we will have a chiral 
spin liquid with a weak ferromagnetic moment. Figure 2 
gives two schematic phase diagrams in the x — T plane. 

Important questions as to why ferromagnetism occur 
in a narrow range of doping x, the complex transport, 
and sharpening our various estimates and heuristic argu- 
ments by many body methods remains to be done. 

I thank P.W. Anderson for an encouraging discussion. 
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